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INTRODUCTION 


The  Isothermal  elastic-plastic  problems  of  thick-walled  cylinders 
subjected  to  mechanical  and/or  thermal  loadings  have  been  solved  by  many 
Investigators  based  on  different  theories  or  methods. The  yield  stress  In 
all  Isothermal  theories  Is  assumed  to  be  temperature-independent.  Although 
good  progress  has  been  made  recently  in  developing  constitutive  relations  for 
thermo-elastlc-plastlc  and  time-dependent  inelastic  theories. ^*6  the  research 
effort  in  this  area  has  not  reached  a  state  of  completion.  In  addition,  the 
general  solution  of  thermo-elastlc-plastlc  problems  Is  still  very  difficult 
and  frequently  very  costly. As  a  result,  our  research  has  been  directed 
towards  the  development  of  a  special  purpose  computer  program  for  solving 
thick-walled  cylinder  problems  of  potential  Importance  to  the  Army. 

This  report  shows  a  numerical  approach  for  analyzing  the  thermo-elastlc- 
plastlc  problems  of  thick-walled  cylinders  with  temperature-dependent  yield 
stress.  The  cylinder  Is  subjected  to  a  combination  of  internal  pressure  and 
temperature  variation.  The  material  is  assumed  to  obey  the  von  Mlses*  yield 
criterion,  the  associated  flow  theory,  end  the  Isotropic  hardening  rule.  Some 
numerical  results  for  the  displacements  and  stresses  are  preseated. 

THERM3-ELASTIC-PLASTIC  THEORY 

e 

For  small  displacement  analysis,  the  total  strain-rate  tensor  ejj  Is 
composed  of  corresponding  elastic,  plastic,  and  thermal  components  as  follows: 

cij  “  eij*  +  eijp  +  eijT 


References  are  listed  at  the  end  of  this  report. 
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On  the  basis  of  the  above  assumptions,  we  can  readily  find  that 


.  3  . 

ci.1p  *  “  (eP/o)sjj  (9) 

and 

3  *ij  do  . 

X  -  eP  -  ( - on - TJ/H'  (10) 

2  a  J  3T 

where 

So  u)E 

H'  -  -  - -  ,  aiE  -  3a/3e  (11) 

3eP 

3o/3T  -  3o0/3T  +  eP(3H*/3T)  (12) 


and  oQ  Is  the  initial  yield  stress. 

Substituting  Eqs.  (2)  and  (9)  into  Eq.  (1),  one  obtains  the  general 

e  '  e  e 

constitutive  equations  relating  ctj  to  oy  and  T.  For  numerical  solutions  by 
the  finite-element  method  or  the  finite-difference  method,  it  is  desirable  to 

•  •  a 

find  the  inverse  form  which  relates  ojj  to  cy  and  T.  For  the  isotropic- 

hardening,  thermo-elastic-plastic  theory,  the  explicit  inverse  relationships 

are  given  below  in  a  form  slightly  different  from  that  shown  in  Reference  7. 
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EQUATIONS  FOR  THICK-WALLED  CYLINDER 


For  the  leotropic-hardening,  thermo-elastic-plastic  thick-walled 
cylinders,  the  Incremental  form  of  Eq.  (13)  reduces  to 


where 


Aol  “  dij  Aej  "  Aoi*  ,  i  «  r,9,z 
E  v  1 

dj4  - - ( - f  6j4  -  -  Oj'  04  ' ) 

J  1+v  l-2v  J  s  J 


(15) 


Ea  ( 3a/3T)  °l' 

t - - J  AT 

l-2v  1+H '  /  3G  a 


(16) 


°i'  -  ai  “  om  ,  ora  -  (or+O0+az)/3 

In  the  quast-statlc  state  with  no  body  forces,  the  radial  and  tangential 
stresses  must  satisfy  the  equilibrium  equation, 

r(3or/3r)  «  oq  -  or  (17) 

and  the  corresponding  strains  must  satisfy  the  compatibility  equation 

r(3ce/3r)  »  er  -  e©  (18) 

Consider  a  thick-walled  cylinder  of  inner  radius  a  and  external  radius  b.  The 
cylinder  Is  subjected  to  Inner  pressure  and  temperature  (p  and  Ta) ,  external 
pressure  and  temperature  (q  and  Tb) ,  and  end  force  (f).  The  boundary 
conditions  for  the  generalized  plane-strain  conditions  are 


or(a,t)  -  -p  , 

T(a,t)  -  Ta 

(19) 

or(b,t)  -  -q  , 

T(b,t)  -  Tb 

(20) 

,b 

2irJ  rojdr  - 

a 

w*a2p  +  f 

(21) 

where  m  is  0  or  1  for  open-end  or  closed- 

-end  conditions,  respectively. 

The 

temperature  distribution  must  satisfy  the  heat  conduction  equation  subjected 
to  boundary  conditions  (19)  and  (20), 


f 


4 


The  sue  of  elastic  and  thermal  atraln-rataa  la  assumed  to  ba  dataralnad  by  the 
Duhsmal-Haunann  law, 

*lj*  +  eljT  "  B“lI(l+v)  Ofck)  +  aT«ij  (2) 

In  which  B  la  Young's  modulus,  v  is  Poisson's  ratio,  a  is  tha  thermal 
expansion  coefficient,  T  la  tha  rata  of  temperature  change,  is  tha 
Kronecker  delta,  and  Ojj  Is  tha  stress  tensor. 

Tha  plastic  strain-rate  ejjP  Is  derivable  from  tha  plastic  potential 
g(oij)  by  tha  normality  condition 

eijp  -  X  3g/3oij  (3) 

a 

where  X  Is  a  positive  scalar  varlabla. 

The  yield  function  for  non-lsothermal  Isotropic  strain-hardening  material 
can  ba  written  as 


where 


F  -  f(otJ)  -  o  <eP,T) 

eP  -  /  eP  dt 
eP  -  (|  «ijP  i1jP)l/2 


(4) 

(5) 

(6) 


and  o(eP,T)  represents  the  dependence  of  yield  stress  on  the  accumulated 
Increments  of  effective  plastic-strain  and  temperature.  Whan  the  von  Mlses' 
yield  criterion  and  associated  flow  rule  are  adopted, 

f  "  *  "  <1  *1J  »1J>1/2  <7> 

end  ‘ 


1 

•1J  *  °1J  “  ”  °kk6lj 


(8) 
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where  k  denotes  thermal  dlffusivlty.  For  the  special  case  of  steady  state 
distribution,  the  temperature  is  given  by 

T  -  Ta  +  (Tb-Ta>  log(r/a)/log(b/a)  (23) 

INCREMENTAL  FINITE-DIFFERENCE  FORMULATIONS 

For  loading  beyond  the  elastic  limit,  an  Incremental  approach  of  the 
finite-difference  formulation  is  used.  The  cross-section  of  the  tube  is 
divided  into  n  rings  with  ri*a,r2,..*»  r^-p, . . . ,rn-n“b,  where  p  is  the  radius 
of  the  elastic-plastic  Interface.  At  the  beginning  of  each  increment  of 
loading,  the  distribution  of  temperature,  displacements,  strains,  and  stresses 
is  assumed  to  be  known  and  we  want  to  determine  Au,  Aer,  Acq,  Aez,  A or,  Aog, 
Aoz  at  all  grid  points  for  the  applied  Incremental  loading,  Ap,  Aq,  Af,  ATi 
(1  ■  1  to  nH).  Since  the  Incremental  stresses  are  related  to  the  Incremental 
strains  by  the  Incremental  form  (Eq.  (13))  and  Au  -  rAeg,  there  exist  only 
three  unknowns  at  each  station  that  have  to  be  determined  for  each  Increment 
of  loading.  Accounting  for  the  fact  that  the  axial  strain  ez  is  Independent 
of  r,  the  unknown  variables  in  the  present  formulation  are  (Aco)^,  (Aer)^,  for 
1  -  1 ,2,. . .n,n+l ,  and  Aez. 

The  equation  of  equilibrium  (17)  and  the  equation  of  compatibility  (18) 
are  valid  for  both  the  elastic  and  the  plastic  regions  of  a  thick-walled  tube. 
The  finite-difference  forms  of  these  two  equations  at  1  -  l,...,n  are  given  by 

ci(Aor)i  +  C2(Aoe)i  +  c3(Aor)i+i  -  C5  (24) 

and 

ci(Aee)i  +  C2(Aer)i  +  C3(Acq)i+i  -  04  (25) 
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where 


C1  •  ri+l  -  2ri  »  c2  “  -ri-K  ♦  *i  »  *3  •  *1 
c«  ■  (ri+l-ri)(er-ee)i  -  ri[(ee)i+i  -  (ee>il 

c5  “  (*i+i"ri)(°0“°r)i  ”  riK°r)i+l  ~  (<>r>il  (26) 

Substitution  of  tbs  incremental  stress-strain  relations  (IS)  into  Eq.  (24) 
leads  to 

*6(669)1  +  C7(Aer)i  +  C8(Ac0)i+i  +  C9(Aer)1+1  +  cio  Aeg  -  cn  (27) 

where 

eg  -  ci(di2)l  +  C2(d22>i  *  *8  “  *3(d12>i+l 
*7  ■  *l(dll)i  +  *2(d21>i  .  *9  *  *3(dll)i+l 
cio  ■  *l<d13>i  +  *2(d23)i  +  *3(d13>i+l 
cn  -  ci(Aor*)i  +  C2(Aoe*)i  +  C3(Aor*)iTi 

+  C2(or-oe)i  +  *3l(°r)i  "  (°r)i+ll  (28) 

The  finite-difference  forms  of  the  boundary  conditions  (19),  (20),  and  (21) 
are 

(di2)l(6ee)i  +  (dll)l(Acr)i  +  (d 1 3)  1  Aee  -  -Ap  +  (Aor*)i  (29) 

(d12)tt+l(6£e)n4-l  +  (dll)ttfl(6er)n+i  +  <d13>iH-l  Aez  m  -M  +  (Aor*>n+l  (30) 

and 

n 

1  [ci21(Ae9)i  +  ci3i(Aer)i  +  ci41(AE0)i+1  +  ci51(Aer)i+1] 

1-1 

n  n 

+  (  [  cie1^6*  *  wa2Ap  +  Af/w  +  l  C171  (31) 

i-1  1-1 
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EQUATIONS  FOR  THICK-WALLED  CYLINDER 


For  the  leotropic-hardening,  thermo-elastic-plastic  thick-walled 
cylinders,  the  Incremental  form  of  Eq.  (13)  reduces  to 


where 


Aol  “  dij  Aej  "  Aoi*  ,  i  «  r,9,z 
E  v  1 

dj4  - - ( - f  6j4  -  -  Oj'  04  ' ) 

J  1+v  l-2v  J  s  J 


(15) 


Ea  ( 3a/3T)  °l' 

t - - J  AT 

l-2v  1+H '  /  3G  a 


(16) 


°i'  -  ai  “  om  ,  ora  -  (or+O0+az)/3 

In  the  quast-statlc  state  with  no  body  forces,  the  radial  and  tangential 
stresses  must  satisfy  the  equilibrium  equation, 

r(3or/3r)  «  oq  -  or  (17) 

and  the  corresponding  strains  must  satisfy  the  compatibility  equation 

r(3ce/3r)  »  er  -  e©  (18) 

Consider  a  thick-walled  cylinder  of  inner  radius  a  and  external  radius  b.  The 
cylinder  Is  subjected  to  Inner  pressure  and  temperature  (p  and  Ta) ,  external 
pressure  and  temperature  (q  and  Tb) ,  and  end  force  (f).  The  boundary 
conditions  for  the  generalized  plane-strain  conditions  are 


or(a,t)  -  -p  , 

T(a,t)  -  Ta 

(19) 

or(b,t)  -  -q  , 

T(b,t)  -  Tb 

(20) 

,b 

2irJ  rojdr  - 

a 

w*a2p  +  f 

(21) 

where  m  is  0  or  1  for  open-end  or  closed- 

-end  conditions,  respectively. 

The 

temperature  distribution  must  satisfy  the  heat  conduction  equation  subjected 
to  boundary  conditions  (19)  and  (20), 


f 
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where 


Ari  ■  ri+1  -  *1  .  C121  “  (Ari)ti(d23>i 

C131  -  (Ar1)ri(di3)i  ,  C141  -  (Ari+i)ri+i(d23)i+l 

c15i  *  (Ari)ri^i(d i3>i 

C161  ■  (Ari){ri(d33)i  +  r1+i(d33>1+1} 

cn1  -  (Ari)(ri(Aoz#)i  +  ri+i( Aos*)i+i)  (32) 

Now  we  can  form  a  system  of  2n+3  equations  for  solving  2n+3  unknowns,  (Aeg)i, 

(Aer)it  at  i  -  1,2 . n,n+l  and  Aez.  Equations  (29),  (30),  and  (31)  are 

taken  as  first  and  last  two  equations,  respectively,  and  the  other  2n 
equations  are  set  up  at  1  -  l,2,...,n  using  equations  (25)  and  (27).  The 
final  system  is  an  unsymmetrlc  matrix  of  arrow  type  with  the  nonzero  terms 
appearing  in  the  last  row  and  column  and  others  clustered  about  the  main 
diagonal,  two  below  and  one  above. 

NUMERICAL  RESULTS  AND  DISCUSSIONS 

The  numerical  results  for  two  particular  problems  follow.  The  first 
problem  is  a  closed-end  thick-walled  cylinder  subjected  to  varying  Internal 
pressure  p  and  temperature  T  as  shown  in  Figure  1.  The  heating  is  uniform 
throughout  the  thickness,  but  the  initial  yield  stress  is  temperature- 
dependent  as  shown  in  Figure  1.  The  other  material  constants  are  E  ■  86,666 
psi,  v  -  0.3,  u  *  0.0,  a  -  0.0.  The  numerical  results  for  the  radial 
displacements  (ua  and  %)  at  the  inside  and  outside  surface  are  shown  in 
Figure  2.  The  percentage  of  plastic  sons  is  also  shown  in  the  figure  by  the 
dotted  line.  The  entire  cylinder  is  elastic  during  the  time  Interval  10  to 
12.  The  results  for  the  three  stress  components  at  selected  time  t  ■  4,  8, 
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where  k  denotes  thermal  dlffusivlty.  For  the  special  case  of  steady  state 
distribution,  the  temperature  is  given  by 

T  -  Ta  +  (Tb-Ta>  log(r/a)/log(b/a)  (23) 

INCREMENTAL  FINITE-DIFFERENCE  FORMULATIONS 

For  loading  beyond  the  elastic  limit,  an  Incremental  approach  of  the 
finite-difference  formulation  is  used.  The  cross-section  of  the  tube  is 
divided  into  n  rings  with  ri*a,r2,..*»  r^-p, . . . ,rn-n“b,  where  p  is  the  radius 
of  the  elastic-plastic  Interface.  At  the  beginning  of  each  increment  of 
loading,  the  distribution  of  temperature,  displacements,  strains,  and  stresses 
is  assumed  to  be  known  and  we  want  to  determine  Au,  Aer,  Acq,  Aez,  A or,  Aog, 
Aoz  at  all  grid  points  for  the  applied  Incremental  loading,  Ap,  Aq,  Af,  ATi 
(1  ■  1  to  nH).  Since  the  Incremental  stresses  are  related  to  the  Incremental 
strains  by  the  Incremental  form  (Eq.  (13))  and  Au  -  rAeg,  there  exist  only 
three  unknowns  at  each  station  that  have  to  be  determined  for  each  Increment 
of  loading.  Accounting  for  the  fact  that  the  axial  strain  ez  is  Independent 
of  r,  the  unknown  variables  in  the  present  formulation  are  (Aco)^,  (Aer)^,  for 
1  -  1 ,2,. . .n,n+l ,  and  Aez. 

The  equation  of  equilibrium  (17)  and  the  equation  of  compatibility  (18) 
are  valid  for  both  the  elastic  and  the  plastic  regions  of  a  thick-walled  tube. 
The  finite-difference  forms  of  these  two  equations  at  1  -  l,...,n  are  given  by 

ci(Aor)i  +  C2(Aoe)i  +  c3(Aor)i+i  -  C5  (24) 

and 

ci(Aee)i  +  C2(Aer)i  +  C3(Acq)i+i  -  04  (25) 
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10,  and  13  ara  shorn  In  Figures  3  through  3.  The  differences  for  the 
displacements  end  stresses  et  t  -  4  end  8  cleerly  demonstrate  the  effect  of 


temperature-dependence  of  the  yield  stress.  The  seee  problem  with  the 
plane-strain  condition  has  been  solved  by  the  Automatic  Dynamic  Incremental 
Analysis  (ADINA)  program.8  For  comparison  purposes,  the  AD  I  HA  results  for  the 
radial  displacement  at  the  outside  and  the  residual  stress  distribution 
through  the  wall  at  time  t  -  10  are  elso  shown  in  Figures  2  and  4.  The  agree¬ 
ment  is  excellent  for  the  stresses  and  good  for  the  displacement.  The  small 
differences  in  the  displacement  response  may  be  due  to  the  end  conditions  and 
the  methods  of  approaches.  The  numerical  results  reported  here  ere  based  on 
the  finite-difference  formulations  with  n  -  100. 

As  a  second  example,  let  us  consider  a  closed-end  tube  subjected  to  inner 
temperature  Ta  only.  The  numerical  results  mere  based  on  the  following 
parameters:  b  -  2",  a  -  1",  n  -  100,  v  •  0.3,  B  ■  30  x  10s  psi,  w  •  0.0,  o  - 
7.75  x  10-6  ln./ln./*F,  oQ  -  30  x  103  psi,  o/oq  -  1.0  -*/2  x  1<T3/#F.  When 
the  temperature  gradient  is  of  sufficient  magnitude,  yielding  will  first 
expand  from  the  inside.  At  larger  temperature  gradient,  the  plastic  zone  will 
expand  from  both  the  Inside  and  outside  surface  toward  the  interior.  The 
relation  between  the  inside  temperature  and  elastic-plastic  Interface  is  shown 
in  Figure  6.  The  stresses  in  e  closed-end  cylinder  subjected  to  temperature 
gradient  of  400*F  are  shown  in  Figure  7.  The  special  case  when  the  yield 
stress  is  assumed  to  be  temperature-independent  was  considered  in  an  earlier 
paper. *  Those  earlier  results  ere  shown  in  Figures  6  and  7  by  the  dotted 
lines.  A  comparison  of  the  results  between  the  solid  and  dotted  lines  shows 
the  effects  of  temperature-dependence  of  the  yield  stress. 
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where 


C1  •  ri+l  -  2ri  »  c2  “  -ri-K  ♦  *i  »  *3  •  *1 
c«  ■  (ri+l-ri)(er-ee)i  -  ri[(ee)i+i  -  (ee>il 

c5  “  (*i+i"ri)(°0“°r)i  ”  riK°r)i+l  ~  (<>r>il  (26) 

Substitution  of  tbs  incremental  stress-strain  relations  (IS)  into  Eq.  (24) 
leads  to 

*6(669)1  +  C7(Aer)i  +  C8(Ac0)i+i  +  C9(Aer)1+1  +  cio  Aeg  -  cn  (27) 

where 

eg  -  ci(di2)l  +  C2(d22>i  *  *8  “  *3(d12>i+l 
*7  ■  *l(dll)i  +  *2(d21>i  .  *9  *  *3(dll)i+l 
cio  ■  *l<d13>i  +  *2(d23)i  +  *3(d13>i+l 
cn  -  ci(Aor*)i  +  C2(Aoe*)i  +  C3(Aor*)iTi 

+  C2(or-oe)i  +  *3l(°r)i  "  (°r)i+ll  (28) 

The  finite-difference  forms  of  the  boundary  conditions  (19),  (20),  and  (21) 
are 

(di2)l(6ee)i  +  (dll)l(Acr)i  +  (d 1 3)  1  Aee  -  -Ap  +  (Aor*)i  (29) 

(d12)tt+l(6£e)n4-l  +  (dll)ttfl(6er)n+i  +  <d13>iH-l  Aez  m  -M  +  (Aor*>n+l  (30) 

and 

n 

1  [ci21(Ae9)i  +  ci3i(Aer)i  +  ci41(AE0)i+1  +  ci51(Aer)i+1] 

1-1 

n  n 

+  (  [  cie1^6*  *  wa2Ap  +  Af/w  +  l  C171  (31) 

i-1  1-1 
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This  report  presents  the  numerical  results  for  a  closed-end  cylinder 
subjected  to  varying  internal  pressure  and/or  temperature.  The  thermal 
problem  is  due  to  uniform  heating  or  thermal  gradient.  If  the  temperature 
distribution  is  solved  by  a  transit  analysis,  the  corresponding  thermal 
stresses  can  be  calculated.  The  result  of  this  transit  thermal  problem  is  to 
be  reported  at  a  later  date. 
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throughout  the  thickness,  but  the  initial  yield  stress  is  temperature- 
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Figure  6.  Blastic-Plaatic  Interface  vs.  Temperature 
Gradient. 
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